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The  Accelerator  Driven  Subcritical  Reactor  (ADSR)  is  one  of  the  reactor  designs  proposed  for  future 
nuclear  energy  production.  Interest  in  the  ADSR  arises  from  its  enhanced  and  intrinsic  safety  character¬ 
istics,  as  well  as  its  potential  ability  to  utilize  the  large  global  reserves  of  thorium  and  to  burn  legacy  acti¬ 
nide  waste  from  other  reactors  and  decommissioned  nuclear  weapons.  The  ADSR  concept  is  based  on  the 
coupling  of  a  particle  accelerator  and  a  subcritical  core  by  means  of  a  neutron  spallation  target  interface. 
One  of  the  candidate  accelerator  technologies  receiving  increasing  attention,  the  Fixed  Field  Alternating 
Gradient  (FFAG)  accelerator,  generates  a  pulsed  proton  beam.  This  paper  investigates  the  impact  of 
pulsed  proton  beam  operation  on  the  mechanical  integrity  of  the  fuel  pin  cladding.  A  pulsed  beam 
induces  repetitive  temperature  changes  in  the  reactor  core  which  lead  to  cyclic  thermal  stresses  in  the 
cladding.  To  perform  the  thermal  analysis  aspects  of  this  study  a  code  that  couples  the  neutron  kinetics 
of  a  subcritical  core  to  a  cylindrical  geometry  heat  transfer  model  was  developed.  This  code,  named  PTS- 
ADS,  enables  temperature  variations  in  the  cladding  to  be  calculated.  These  results  are  then  used  to  per¬ 
form  thermal  fatigue  analysis  and  to  predict  the  stress-life  behaviour  of  the  cladding. 

©  2011  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Accelerator  Driven  Subcritical  Reactors  (ADSR)  are  one  of  the 
possible  future  fission  reactor  systems  under  consideration  that 
target  sustainability,  enhanced  safety  and  economic  competitive¬ 
ness.  In  addition,  ADSRs  can  be  used  to  incinerate  problematic  acti¬ 
nide  waste  and  long-live  radiotoxic  fission  products.  Interest  in 
ADSRs  also  arises  from  the  possibility  of  deploying  thorium  in 
them.  The  use  of  thorium  in  an  ADSR  was  first  proposed  by  Bow¬ 
man  (1998).  Thorium  is  three  to  four  times  more  abundant  in  nat¬ 
ure  than  uranium  and  under  certain  circumstances  produces  fewer 
minor  actinides  (International  Atomic  Energy  Agency,  2005). 

The  concept  of  an  ADSR  is  based  on  the  coupling  of  a  particle 
accelerator  that  delivers  a  beam  of  protons  of  about  1  GeV  energy 
to  initiate  spallation  reactions  in  a  heavy  metal  target  and  a  sub¬ 
critical  reactor  fed  by  the  spallation  neutrons.  The  fact  that  the 
reactor  core  is  inherently  subcritical  means  that  the  reactor  can 
be  shut  down  simply  by  shutting  off  the  accelerator.  Among  the 
accelerator  technologies  under  consideration  for  deployment  in 
an  ADSR,  the  Fixed  Field  Alternating  Gradient  (FFAG)  concept 
seems  to  be  a  promising  candidate  (Bungau  et  al.,  2008).  The  FFAG 
offers  better  performance  than  a  cyclotron  in  terms  of  its  beam 
characteristics  such  as  beam  focusing  and  rapid  acceleration  (Mori, 
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2006).  Moreover,  FFAGs  have  a  lower  construction  cost  than  linear 
accelerators  (Takahashi,  2000).  The  FFAG  was  developed  for  rapid 
acceleration,  which  is  ideal  for  systems  requiring  both  high  beam 
power  and  high  repetition  rates  (Symon  et  al.,  1956).  An  FFAG 
accelerator  has  already  been  coupled  to  a  reactor  core  in  Kyoto 
University  Research  Reactor  Institute  (KURRI)  and  studies  of  its 
performance  conducted  (Yamamoto  and  Shiroya,  2003;  Tanigaki 
et  al.,  2006).  An  FFAG  delivers  accelerated  protons  in  the  form  of 
discontinuous  bunches,  which  in  an  ADSR  would  result  in  a  pulsed 
spallation  source.  This,  in  turn,  would  lead  to  oscillations  in  the 
number  of  fission  reactions  within  the  subcritical  core  and  there¬ 
fore  in  the  core  power.  Temperatures  in  the  core  can,  in  conse¬ 
quence,  also  be  expected  to  oscillate. 

This  paper  presents  work  done  in  coupling  a  fuel  pin  heat  trans¬ 
fer  model  to  a  neutron  kinetics  model  in  order  to  study  the  dy¬ 
namic  response  of  the  fuel  to  pulsed-beam  operation  in  an  ADSR. 
This  coupling  was  required  to  facilitate  thermo-mechanical  stress 
analyses,  based  on  which  the  consequences  of  various  accelerator 
beam  transients  can  be  quantified.  The  pulsed-beam  delivery  in 
FFAG  accelerators  will  induce  cyclic  thermal  stresses  in  the  fuel 
pellets  and  cladding.  This  study  focusses  attention  on  the  ther¬ 
mo-mechanical  response  of  the  fuel  cladding  as  this  is  the  barrier 
that  prevents  fission  fragments  from  leaking  into  the  coolant 
circuit.  As  long  as  the  cladding  maintains  its  integrity  thermo¬ 
mechanical  damage  to  the  fuel  pellets  is  of  limited  significance. 
The  thermally  induced  stress  variation  in  the  cladding  will 
determine  the  number  of  allowable  cycles,  both  in  normal 
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Nomenclature 

P 

coolant  density  (kg  m-3) 

KP 

fatigue  plasticity  coefficient 

Pc 

clad  density  (kg  m-3) 

m 

coolant  mass  flow  rate  (kg  s_1) 

pf 

fuel  density  (kg  m-3) 

Hc 

number  of  radial  clad  nodes 

Oa 

stress  amplitude 

nf 

number  of  radial  fuel  nodes 

y 

A 

fatigue  strength  coefficient 

q' 

linear  heat  flux  (W  m_1) 

channel  flow  area  (m2) 

q" 

surface  heat  flux  (W  m-2) 

b 

Basquin’s  exponent 

q'" 

volume  heat  flux  (W  m-3) 

cP 

coolant  specific  heat  capacity  (J  kg-1 1<-1) 

rc 

clad  outer  radius  (m) 

Cpc 

clad  specific  heat  capacity  (J  kg-1  K1) 

rf 

fuel  pellet  outer  radius  (m) 

CPf 

fuel  specific  heat  capacity  (J  kg-1  K-1) 

rg 

clad  inner  radius  (m) 

hw 

clad-coolant  heat  transfer  coefficient  (W  m-2  K-1) 

Tb 

coolant  bulk  temperature  (K) 

hg 

gap  heat  transfer  coefficient  (W  m-2  I<-1) 

Tc 

clad  outer  surface  temperature  (K) 

K 

cyclic  strength  coefficient 

kc 

clad  thermal  conductivity  (W  K-1  m-1) 

kf 

fuel  thermal  conductivity  (W  I<-1  m_1) 

pulsed-beam  operation  and  in  the  face  of  the  larger  transients  in¬ 
duced  by  beam  trips.  The  coupling  between  models  of  the  reactor 
neutron  kinetics  and  fuel  pin  heat  transfer  was  achieved  through 
the  development  of  the  PTS-ADS  code.  This  code  is  described  in 
the  next  section. 

2.  The  PTS-ADS  code 

PTS-ADS  is  a  Fortran77  code  developed  to  couple  the  neutron 
kinetics  of  subcritical  systems  to  a  fuel  pin  heat  transfer  model. 
The  code  uses  an  implicit  finite  difference  method,  based  on  for¬ 
ward  time-marching.  At  present,  the  code  models  a  single  channel 
in  a  fast  neutron  spectrum  reactor.  The  PTS-ADS  reactor  kinetics 
are  described  by  a  six-group  point  kinetics  model  with  tempera- 
ture-to-reactivity  feedback  included.  The  single  channel  heat 
transfer  model  consists  of  a  fuel  heating  model  and  a  heat  transfer 
model  from  a  heated  wall  to  the  coolant  bulk. 

The  reactor  physics  of  an  ADSR  is  fundamentally  different  from 
that  of  critical  reactors.  Its  subcriticality  has  many  implications  for 
the  neutronic  behaviour  of  the  core  and  thus  for  reactor  control. 
While  the  power  level  of  critical  reactors  depends  on  the  fuel  mass 
present  in  the  core  and  the  deployment  of  control  materials  (control 
rods,  soluble  boron  etc.),  the  power  level  in  an  ADSR  depends  on  the 
margin  of  subcriticality  and  the  accelerator  beam  current.  In  subcrit¬ 
ical  reactors  the  total  neutron  flux  is  the  summation  of  the  source 
neutrons,  prompt  fission  neutrons  and  delayed  neutrons.  The  source 
neutrons  are  produced  via  the  spallation  process  through  the  inter¬ 
action  between  the  high  energy  protons  in  the  accelerator  beam  and 
the  heavy  metal  spallation  target,  and  therefore  their  population  is 
independent  of  the  multiplying  medium  characteristics.  In  contrast, 
the  populations  of  prompt  and  delayed  neutrons  depend  of  the  prop¬ 
erties  of  the  multiplying  medium  and  the  intensity  of  the  source.  The 
neutronic  kinetics  can  be  described  by  a  six-group  point  kinetics 
model  with  a  time-dependent  neutron  source.  The  use  of  a  point 
kinetics  model  to  describe  the  physics  of  ADSRs  was  successfully 
demonstrated  by  Eriksson  et  al.  (2005). 

The  nominal  intensity  of  the  source  term  when  the  accelerator 
beam  is  on  can  be  calculated  based  on  steady-state  operating  con¬ 
ditions  and  power  normalisation.  When  the  beam  is  off,  the  rate  of 
decay  of  the  neutron  population  depends  strongly  on  the  type  of 
coolant,  i.e.  whether  the  ADSR  is  operating  with  a  fast  or  thermal 
neutron  spectrum.  This  translates  into  different  neutron  lifetimes 
in  the  point  kinetics  equations. 

In  general,  temperature-to-reactivity  feedback  effects  are  not  as 
important  in  an  ADSR  as  in  conventional  critical  reactors  due  to 
the  presence  of  the  source  neutrons:  these  reduce  the  ADSR’s  sensi- 


Fig.  1.  The  PTS-ADS  thermo-hydraulic-neutronic  coupled  model. 

tivity  to  reactivity  changes.  However,  when  the  system  is  close  to 
criticality,  reactivity  feedback  mechanisms  become  more  signifi¬ 
cant.  Two  such  mechanisms  have  been  included  in  the  PTS-ADS  neu¬ 
tronic  model:  fuel  temperature  feedback  (Doppler  feedback)  and 
coolant  density  feedback.  The  former  depends  on  the  fuel  design 
and  composition;  the  latter  depends  on  the  choice  of  coolant  and 
flow  parameters. 

2A.  PTS-ADS  physical  modelling 

A  flow  chart  of  the  PTS-ADS  coupled  model  is  shown  in  Fig.  1.  A 
run  starts  by  calculating  the  power  generated  in  the  fuel  from  the 
initial  conditions  and  the  temperature  distribution  in  the  fuel  and 
clad  is  obtained  by  applying  the  heat  conduction  model.  The  corre¬ 
sponding  Doppler  reactivity  feedback  is  calculated  and  saved.  The 
coolant  bulk  temperature  is  then  calculated  using  the  convective 
heat  transfer  model  with  the  appropriate  heat  transfer  coefficient 
correlation.  Having  determined  the  coolant  temperature,  the  cool¬ 
ant  temperature-to-reactivity  feedback  can  then  be  calculated. 
Both  reactivity  feedbacks  are  then  inserted  into  the  point  kinetics 
model  which  calculates  the  power  generation  in  the  fuel,  and  so  on. 

2AA.  Fuel  pin  heat  transfer  model 

The  PTS-ADS  heat  transfer  model  uses  the  following  assump¬ 
tions  in  the  computation  of  the  temperature  distribution: 

•  The  axial  power  distribution  has  a  sinusoidal  form.  1 


1  The  power  distribution  in  an  ADSR  depends  on  the  core’s  effective  multiplication 
factor,  but  the  distribution  can  be  assumed  to  follow  the  fundamental  harmonic  if  the 

system  is  not  highly  subcritical  (Kadi  and  Carminati,  2002). 
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Fig.  2.  Fuel  rod  geometry. 

•  The  heat  generation  rate  is  uniform  radially  across  the  fuel  pin. 

•  Heat  transfer  only  occurs  in  the  radial  direction. 


The  fuel  pin  cylindrical  geometry  is  shown  in  Fig.  2.  At  each  ax¬ 
ial  step,  the  model  computes  the  radial  temperature  distribution 
by  solving  the  heat  transfer  equation,  which,  in  cylindrical  coordi¬ 
nates,  can  be  written  as: 


pcP 


&r 

dt 


\d_ 

rdr 


(1) 


where  Q.  is  the  heat  source  term,  k  is  the  thermal  conductivity,  Cp  is 
the  specific  heat  capacity  and  p  is  the  density  of  the  material. 

The  boundary  conditions  are  given  by: 


(“kcf)r  r  =MTc(0-™) 


(2) 

(3) 

(4) 

(5) 


where  the  various  parameters  are  defined  in  the  Nomenclature. 

To  solve  the  radial  heat  transfer  equation  subject  to  the  defined 
set  of  boundary  conditions  (Eqs.  (2)-(5))  the  finite  difference  meth¬ 
od  is  used  to  discretise  both  the  time  and  space  derivatives.  The 
space  discretisation  is  achieved  by  constructing  a  radial  mesh  as 
shown  in  Fig.  3. 

Using  standard  methods  Eq.  (1)  can  be  discretised  giving  the 
coefficients  for  a  standard  matrix  representation  as  detailed  in 
Appendix  A. 


b\  Ci 

a2  b2  c2 


A  = 


Oi  bi  ^ 


Q)if+nc- 1 


bns+nc-\ 

Q-nf+ric 


Crif+nc-1 

brif+ric 


-  T,  ' 

-  -q'"/(pfCPf)  - 

Tn/ 

-q'"/(PfCpf) 

Trif+ 1 

;  B  = 

0 

Trif+ric- 1 

0 

Tnf+ric 

_  ~hwTb/ (pcCPc)_ 

(8) 


(9) 


where  the  various  parameters  are  defined  in  the  Nomenclature.  The 
coefficients  of  matrix  A  are  detailed  in  Appendix  A. 

Since  the  coefficients  of  matrix  A  are  variable,  Eq.  (7)  is  solved 
using  Newton’s  iteration  method.  The  solution  starts  with  an  initial 
guess  of  the  radial  temperature  distribution  (for  instance,  it  can  be 
assumed  to  be  uniform  and  equal  to  the  coolant  inlet  temperature), 
then  the  second  distribution  is  calculated  based  on  the  first,  the 
third  based  on  the  second  and  so  on.  The  calculation  is  halted  when 
convergence  is  achieved.  At  each  iteration,  Eq.  (7)  is  solved  using 
the  tridiagonal  matrix  algorithm,  a  simplified  form  of  Gaussian 
elimination  (Conte  and  de  Boor,  1980). 


2.1.3.  Transient  numerical  solution 

Because  of  its  better  numerical  stability,  the  implicit  method  is 
used  for  the  transient  case.  The  time  term  in  differenced  form  at 
node  i  is: 


vr 

dt 


TJ+1  _  TJ 
1  i _ 21 

At 


+  O(At) 


The  matrix  form  of  Eq.  (1)  can  be  written  as: 


(10) 


P+l  =  P  +  A  t(4J+1P+1  +  g+1 ) 


(11) 


Rearranging,  Eq.  (11)  can  be  expressed  as: 


2.1.2.  Steady-state  numerical  solution 

In  steady-state  operation,  Eq.  (1)  reduces  to: 


1  1  d 

pCp  r  dr 


OT 


1 


r\C~—  +  — =  0 


dr)  pCp 


4'P+I  =  B' 


where 


(6) 


4'  =  /-  At4,+1 


Using  the  finite  difference  forms,  Eq.  (6)  becomes  in  matrix-vector 
form: 


B'  =  P  +  AfB,+1 


(12) 


(13) 

(14) 


41  =  B  (7) 

where  A  is  a  tridiagonal  matrix  of  temperature-dependent  coeffi¬ 
cients,  T  is  the  vector  of  radial  temperatures  and  B  is  the  radial  heat 
generation  vector: 


■  gap  ■ 

d. 


cladding 

<4 


coolant 


1  2  3 


itf  rif  + 1  rij-  +  2 


with  /  being  the  identity  matrix. 

A  and  B  are  the  same  as  in  the  steady-state  case.  Eq.  (12)  is  of 
the  same  form  as  Eq.  (7)  and  can  therefore  be  solved  in  a  similar 
manner. 

The  coolant  bulk  temperature  Tb  is  found  by  solving  the  follow¬ 
ing  energy  equation: 

c'(',4§+m§)-,'<2>  1,51 

where  the  various  parameters  are  defined  in  the  Nomenclature, 
with  q'(z),  the  linear  heat  flux,  being: 

q'(z)  =  2nrchw(Tnf+„c  -  Tb)  (16) 


Fig.  3.  Fuel  rod  radial  mesh. 


Eq.  (15)  can  be  written  in  finite  difference  form  as: 
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CppA 


n:/  n., 


Af 


Cpth 


VV 1 


Az 


=  2urchw(vnf+nc  l  -  T #w) 

where  /  is  the  axial  node  index.  Hence,  rearranging  Eq.  (17): 


2nrchw  /  j  j  \ 

rn  f  rpj  r^j 

'll 

_  CppA  \  nr+n‘J  b4 

pAAz\b’1 

vj 

(17) 


(18) 


2.1.4.  Neutronics  model 

PTS-ADS  uses  a  point  kinetics  model  with  6  groups  of  delayed 
neutrons  to  calculate  the  variation  of  power  with  time  during  tran¬ 
sients.  The  mathematical  formulation  of  the  model  is: 


om 

dt 


+S(t) 

71  i=l 


(19) 


^  =  (20) 

where  P(t)  is  the  normalised  prompt  power,  Q(t)  is  the  population 
of  precursor  group  i ,  p  is  the  effective  delayed  neutron  fraction,  ft 
is  the  delayed  neutron  fraction  of  group  i ,  p  is  the  total  reactivity, 
A  is  the  prompt  neutron  generation  time,  and  S(t)  is  the  indepen¬ 
dent  source  term. 

The  total  reactivity  can  be  written  as: 

pit)  =  Pi  0)  +  A  PfuetiTf)  +  A  pC00lam(Tb)  (21) 

where  A pfuei  (T/)  is  the  fuel-temperature-to-reactivity  feedback,  also 
known  as  the  Doppler  reactivity  feedback,  and  A pcooiant  {W  is  the 
coolant-temperature-to-reactivity  feedback  due  to  the  variation  in 
coolant  density  with  temperature. 

The  finite  difference  forms  of  Eqs.  (19)  and  (20)  are: 


piA  _  pi 

At 

dp  -  c{ 

Af 


pimipj+i 


+  jfliCp+Si+\t) 

i=  1 


^p,+1  -  2,-cf 1 


which  can  be  rearranged  to  give: 

P'  +AfEti  (Sm)  +At^+1(f) 

f  “l-At^-At2^^) 


(22) 

(23) 


(24) 


Qj+1 


Cj  +  AfAp'+1 
1  +  Atli 


(25) 


Using  Eq.  (24)  the  normalised  power  at  time  step  j  +  1  can  be 
calculated  based  on  parameter  values  from  the  previous  time  step, 
and  then  the  precursor  group  populations  can  be  updated  using  Eq. 
(25). 

The  initial  conditions  for  this  point  kinetics  model  are: 


P1  =  1 

rl=Pft 

'■  Ah 


(26) 

(27) 


ADSR  to  beam  interruptions  predicted  using  a  set  of  ten  different 
codes.  The  three  selected  codes  for  the  comparison  with  PTS-ADS 
are:  TRAC-MOD  (Liles,  1982),  SAS4ADS  (Dagan  et  al.,  2002)  and  SI- 
MADS  (Schikorr,  2002).  These  three  codes  represent  well  the  range 
of  data  scatter  in  the  benchmark  study. 

Fig.  4  compares  the  predictions  by  the  four  codes  of  the  axial 
temperature  distributions  in  the  European  experimental  accelera¬ 
tor-driven  system  (XADS)  and  shows  good  agreement  between 
the  PTS-ADS  steady-state  model  and  the  other  three  codes.  It  is 
clear  that  the  small  discrepancies  between  the  codes  reduce 
greatly  with  distance  from  the  fuel  centreline.  The  focus  in  this 
study  is  on  temperature  variations  in  the  cladding  for  which  all 
the  codes  are  in  good  agreement. 

The  PTS-ADS  neutron  kinetics  model  was  tested  similarly  by 
comparing  its  predictions  of  the  magnitude  of  the  normalised 
power  after  beam  restoration  following  interruptions  of  1,  3,  6 
and  12  s  to  those  in  the  benchmark  study  produced  using  the  three 
selected  codes.  The  comparison  is  shown  in  Fig.  5.  This  shows  that 
compared  to  the  other  three  codes  PTS-ADS  has  slightly  over-esti¬ 
mated  the  peak  normalised  power  reached  after  long  beam  inter¬ 
ruptions.  The  difference,  however,  is  less  than  2%.  Although  we 
have  endeavoured  to  follow  the  benchmark  recommendations  as 
closely  as  possible,  the  discrepancy  in  the  peak  normalised  power 
in  Fig.  5  can  be  attributed  to  the  effects  of  slightly  different  model 
assumptions  or  data.  The  same  cause  of  discrepancies  was  identi¬ 
fied  by  the  participants  of  the  benchmark  exercise  (D’Angelo 
et  al.,  2004). 

The  same  cases  were  used  to  validate  the  PTS-ADS  transient 
thermal  model.  PTS-ADS  predictions  of  transient  temperature  vari¬ 
ations  following  beam  trips  and  restorations  were  compared  to 
those  of  the  three  reference  codes.  Figs.  6  and  7  show  these  com¬ 
parisons  for  the  transient  variation  of  the  peak  fuel  centreline  tem¬ 
perature  and  the  coolant  outlet  temperature  following  a  1  s  beam 
interruption.  In  both  cases  it  can  be  seen  that  the  PTS-ADS  results 
agree  well  with  the  predictions  of  the  reference  codes. 


3.  Thermo -mechanical  response  to  pulsed  operation  transients 


2.2.  PTS-ADS  validation 

The  PTS-ADS  coupled  model  has  been  validated  by  comparing 
its  predictions  for  both  the  steady-state  and  transient  cases  to 
the  predictions  of  three  different  codes  included  in  a  benchmark 
study  on  beam  interruptions  in  an  accelerator-driven  system 
(D’Angelo  et  al.,  2003;  D’Angelo  et  al.,  2004).  The  benchmark  study 
summarises  a  comparative  analysis  of  the  dynamic  response  of  an 


3.1.  Thermal  response 

The  PTS-ADS  code  was  used  to  calculate  the  temperature  varia¬ 
tions  in  fuel  cladding  due  to  pulsed-beam  operation.  The  fuel  pin 
geometry  and  the  physical  properties  of  the  fuel,  cladding  and 
coolant  were  taken  from  the  XADS  fuel  pin  design  (D’Angelo 
et  al.,  2004).  The  values  and  correlations  of  the  thermo-physical 
properties  of  XADS  materials  are  listed  in  Table  1.  Unlike 
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Fig.  6.  XADS  peak  fuel  centreline  temperature  transient  induced  by  a  1  s  beam 
interruption. 
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Fig.  7.  XADS  coolant  outlet  temperature  transient  induced  by  a  1  s  beam 
interruption. 


sodium-cooled  reactors  where  the  temperature  difference  between 
the  cladding  and  the  coolant  bulk  does  not  normally  exceed  10  K, 
the  lower  thermal  conductivity  of  lead  or  lead-bismuth  eutectic 
(LBE),  the  main  coolant  candidates  for  fast  ADSRs,  increases  the 
clad-coolant  temperature  difference  significantly.  Therefore,  an 
appropriate  choice  of  correlation  for  the  heat  transfer  coefficient 
is  very  important.  The  current  version  of  the  PTS-ADS  code  uses 
the  heat  transfer  correlation  derived  by  Mikityuk  (Pfrang  and  Stru- 
we,  2007).  The  adopted  correlation  works  for  triangular  fuel  bun¬ 
dles  and  represents  the  best  fit  of  658  experimental  data  points. 

The  neutronic  dynamic  response  to  pulsed-beam  operation  in 
an  ADSR  depends  on  the  margin  of  subcriticality,  spatial  position 


in  the  core  and  the  beam  characteristics  (Ahmad  and  Parks,  in 
preparation).  For  this  study,  the  effective  multiplication  factor  keff 
is  kept  constant  and  equal  to  0.97276,  the  reference  value  for 
XADS,  for  all  cases  studied.  The  effect  of  the  spatial  position  of 
the  fuel  rod  can  be  represented  by  changing  its  local  power.  The 
accelerator  beam  characteristics,  specifically  the  beam  frequency 
and  beam-off  time,  play  an  important  role  in  shaping  the  neutronic 
and  thermal  dynamic  response  in  the  core.  Since  the  temperature 
variations  in  the  cladding  due  to  pulsed-beam  operation  are  ex¬ 
pected  to  be  small,  the  thermal  cyclic  fatigue  is  expected  to  be  of 
a  high-cycle  nature.  The  influence  of  loading  frequency  is  insignif¬ 
icant  in  the  case  of  high-cycle  fatigue  conditions  in  metals  (Papaky- 
riacou  et  al.,  2001 ).  Consequently,  the  choice  of  the  frequency  is  not 
important  from  a  thermo-mechanical  response  perspective.  All 
cases  were  therefore  run  with  a  frequency  of  1  Hz.  It  should  be 
noted,  however,  that  operating  an  FFAG  accelerator  at  a  frequency 
of  1  Hz  would  mean  a  very  high  instantaneous  peak  current.  One  of 
the  advantages  of  FFAG  accelerators  is  their  ability  to  operate  at 
higher  frequencies  (up  to  1  kHz)  than  conventional  synchrotrons. 
The  use  of  FFAG  technology  for  ADSR  applications  is  a  topic  of  ac¬ 
tive  research  (Tanigaki  et  al.,  2006;  Tanigaki  et  al.,  2007)  and  the 
optimal  values  of  repetition  frequency  and  beam-off  time  have 
not  been  yet  established. 

It  is  expected  that  the  beam-off  time  of  an  FFAG  will  be  of  the 
order  of  10  ms  (Tanigaki  et  al.,  2007).  As  the  beam-off  time  in¬ 
creases,  the  magnitudes  of  the  temperature  fluctuations  in  the 
cladding  inner  and  outer  surface  increase  and  the  effect  of  increas¬ 
ing  the  local  linear  power  becomes  more  significant.  For  this  study 
the  beam-off  time  was  therefore  conservatively  estimated  to  be 
100  ms.  The  linear  power  corresponding  to  that  of  the  midplane  lo¬ 
cal  power  of  the  XADS  fuel  rod  studied  in  (D’Angelo  et  al.,  2003)  is 
91 72  W  m_1.  It  was  estimated  that  for  an  industrial-scale  2  GW(th) 
power  generating  ADSR  the  linear  power  at  midplane  would  be 
25  kW  m_1.  Fig.  8  shows  the  variation  of  the  cladding  inner  surface 
temperature  corresponding  to  this  power  rating  and  a  beam-off 
time  of  100  ms.  It  can  be  seen  that  even  for  this  extreme  case  the 
magnitude  of  the  temperature  oscillation  is  only  about  2  K. 

3.2.  Cladding  thermal  stress 

During  its  service,  the  cladding  may  suffer  from  different  thermo¬ 
mechanical  effects  such  as  creep,  swelling  and  embrittlement  that 
influence  its  mechanical  integrity.  In  the  case  of  a  pulsed-beam 
ADSR,  the  cyclic  variations  of  the  cladding  temperature  will  induce 
thermal  fatigue.  The  significance  of  such  cyclic  behaviour  depends 
on  a  range  of  factors  such  as  the  physical  properties  of  the  cladding 
material,  the  amplitude  and  mean  of  induced  stresses  etc. 

The  total  stress  loading  applied  to  the  cladding  during  operation 
is  the  sum  of  the  mechanical  and  thermal  stresses.  The  mechanical 
stress  oM  can  be  estimated  using  Thin  Cylinder  Theory  (Hearn,  1997) 
as: 


Cm 


(Pgas-Ph)re 

dc 


(28) 


where  pgas  is  the  internal  pressure  due  to  fission  gases,  ph  is  the 
hydrostatic  pressure  applied  on  the  outer  cladding  surface,  rg  is 
the  cladding  inner  radius  and  dc  =  rc-  rg  is  its  thickness.  Assuming 
the  net  pressure  is  1  MPa  (Ceballos-Castillo,  2007),  the  mechanical 
stress  applied  to  the  cladding  is  7.522  MPa.  The  value  of  applied 
mechanical  stress  will  not  greatly  affect  the  thermal  fatigue  calcu¬ 
lations  as  the  stress  variation  is  only  due  to  that  of  the  thermal 
stress. 

The  thermal  stress  acting  on  the  cladding  is  caused  by  the  radial 
temperature  gradient  between  the  inner  and  outer  surfaces.  During 
pulsed  operation,  the  temperature  gradient  will  change  with  time 
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Table  1 

XADS  fuel  pin  material  properties  (D’Angelo  et  al.,  2003). 


Property 

Correlation  or  value 

Clad  thermal  conductivity,  kc 

15.4767  +  0.0034487  W  K1  m”1 

Clad  density,  pc 

7924  kg  rrr3 

Clad  specific  heat  capacity,  Cpc 

620  J  kg”1  K1 

Fuel  thermal  conductivity,  kf 

[0.042  +  2.71  x  1 0— 4T]— 1  +  6.9  x  lO”11!3  W K1  itT1 

Fuel  density, 

10,354  kg  rrr3 

U02  specific  heat  capacity,  CPU02 

81.825  +  0.786957-  1.1552  x  10“372 
+9.9037  x  10_773  -  5.1982  x  10 “^T4 
+1.5241  x  10“137^  -  1.7906  x  10“177^J  kg”1  K1 

Pu02  specific  heat  capacity,  CPPll02 

-4.9236  x  10 6r2  +  240.89  +  0.325567 
-3.5398  x  10_472  +  1.512  x  10“773 
-1.9707  x  10_1174  J  kg-1  K”1 

Fuel  specific  heat  capacity,  CPf 

214.65  r  _i_  55.56  r 

270.21  UPl/02  270.21  U PPu02 

Coolant  density,  p 

11,112  -  1.3757  kg  m“3 

Coolant  specific  heat  capacity,  Cp 

1 46.5  J  kg”1  IC1 

Time(s) 

Fig.  8.  Cladding  inner  surface  temperature  variation  due  to  pulsed  operation  with 
100  ms  beam-off  time  and  1  Hz  frequency  at  25  kW  m_1  linear  power. 


The  stress  amplitude  Acra  =  \  (omax  -  amin)  is  related  to  Nf  by 
Basquin’s  law  (Basquin,  1910): 

Aaa  =  <r}(2  Nf)b  (30) 

where  &f  is  the  fatigue  strength  coefficient  and  b  is  the  fatigue 
strength  exponent  (Basquin’s  exponent).  For  T91  at  300  °C  and 
550  °C  which  define  the  range  of  the  cladding  operating  tempera¬ 
ture,  (T'f  and  b  can  be  obtained  from  the  data  published  in  (Verleene 
et  al.,  2006;  Weisenburger  et  al.,  2008)  which  relate  strain  rather 
than  stress  to  the  number  of  allowable  cycles.  Basquin’s  exponent 
b  would  be  the  same  in  both  forms.  a'f  can  be  written  as  (Tros- 
hchenko  and  Khamaza,  2010): 

a'f  =  K'(Kp)w  (31) 

where  K'  is  the  cyclic  strength  coefficient,  I<p  is  the  fatigue  plasticity 
coefficient,  and  n!  is  the  cyclic  strain  hardening  exponent.  Since  the 
thermal  fatigue  in  the  fuel  cladding  is  of  high-cycle  nature,  the  ef¬ 
fect  of  the  mean  stress  should  be  included  to  improve  accuracy. 
Morrow  proposed  a  modification  to  the  fatigue  strength  coefficient 
to  account  for  the  effect  of  the  mean  stress  (Morrow,  1968): 


and  therefore  induce  compression  and  expansion  stresses.  The 
thermal  stress  gt  applied  to  the  cladding  can  be  estimated  using 
(ASME,  2010): 


(Jt 


EpAT 

2(T^0 


(29) 


Of  =  °f  ~  (32) 

where  am  is  the  mean  stress. 

To  account  for  the  effects  of  local  plasticity  and  creep,  a  correc¬ 
tion  factor  to  the  stress  amplitude  has  to  be  calculated.  The  method 
is  described  in  (ASME,  2010).  The  modified  maximum  equivalent 
strain  range  A emod  can  calculated  as: 


where  E  is  the  Young’s  modulus  of  the  cladding  material,  v  is  its 
Poisson’s  ratio,  p  is  its  coefficient  of  thermal  linear  expansion,  and 
AT  is  the  temperature  difference  between  the  cladding  inner  and 
outer  surfaces.  The  cladding  material  used  in  this  study  is  the  same 
as  that  intended  for  XADS:  9Cr-lMo  (T91)  steel.  T91  steel  is  one  of 
the  proposed  structural  materials  for  ADSR  applications  because  of 
its  good  thermo-mechanical  properties  and  behaviour  under  irradi¬ 
ation  (Rasche  et  al.,  1992).  The  mechanical  properties  of  T91  steel 
used  in  this  study  are  presented  in  Table  2. 

3.3.  Life  prediction  of  fuel  pin  cladding 

The  behaviour  of  fuel  pin  cladding  under  pulsed-beam  opera¬ 
tion  can  be  predicted  by  relating  the  amplitude  of  stress  or  strain 
variations  to  the  number  of  cycles  to  failure  Nf.  Since  the  maximum 
amplitude  of  the  stress  variation  in  all  the  cases  considered  is  low¬ 
er  than  the  yield  stress,  the  fatigue  loads  are  small  and  the  induced 
strain  is  elastic.  Fig.  9  shows  the  variation  in  cladding  hoop  stress 
for  the  same  power  rating  and  beam-off  time  as  used  in  Fig.  8. 


Aemod  =  (ffjK2  Aemax  (33) 

where  K  is  either  the  stress  concentration  factor  or  the  maximum 
value  of  the  theoretical  elastic  stress  concentration  factor,  a*  is 
the  stress  at  the  maximum  strain  range,  g  is  the  stress  at  a  strain 
range  I<Aemax,  and  Aemax  is  the  maximum  equivalent  strain  range. 
For  an  elastic  uniaxial  stress,  I<  can  be  defined  as  the  ratio  of  the 
maximum  stress  to  the  nominal  (mean)  stress,  that  is: 

K=°max  (34) 

(Tm 

The  strain  terms  in  Eq.  (33)  can  be  replaced  by  the  stress  di¬ 
vided  by  Young’s  modulus,  which  is  assumed  to  be  constant  during 
the  stress  cycle  as  the  temperature  change  is  relatively  small.  Eq. 
(33)  can  then  be  written  as: 

Aamod=(^)K2A(Ja  (35) 

and,  as  a  =  I<g *,  simplified  to: 
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Table  2 

T91  steel  mechanical  properties  at  300  °C. 


Property 

Value 

Reference 

Young’s  modulus,  E 

187,900  MPa 

Ashrafi-Nik  (2006) 

Coefficient  of  thermal  linear  expansion,  f 

12.4  x  10-6  K-1 

Ashrafi-Nik  (2006) 

Poisson’s  ratio,  v 

0.3 

Ashrafi-Nik  (2006) 

Yield  stress,  oy 

440  MPa 

Verleene  et  al.  (2006) 

Ultimate  tensile  strength,  aTS 

550  MPa 

Verleene  et  al.  (2006) 

Fig.  9.  Cladding  hoop  stress  variation  due  to  pulsed  operation  with  100  ms  beam- 
off  time  and  1  Hz  frequency  at  25  kW  m-1  linear  power. 


Table  3 

T91  cladding  stress-life  parameters  in  an  LBE-cooled  ADSR. 


Temp  (°C) 

b 

(7 m  (MPa) 

o'f  (MPa) 

<x"  (MPa) 

A  Gmod  (MPa) 

Nf 

300 

-0.052 

29.66 

725.26 

695.60 

0.24 

oo 

550 

-0.093 

69.29 

719a 

649.71 

0.7 

oo 

a  The  fatigue  strength  coefficient  at  550  °C  was  calculated  using  Mitchell  esti¬ 
mation  in  Troshchenko  and  Khamaza  (2010).  The  T91  ultimate  tensile  stress  at 
550  °C  was  taken  from  Chaouadi  (2005). 


A  Gmod  =  KACTa  (36) 

Using  the  temperature  variations  estimated  by  the  PTS-ADS 
code  for  the  defined,  conservative  ADSR  pulsed-beam  operating 
characteristics,  the  stress-life  parameters  of  T91  fuel  cladding  in 
an  LBE  environment  were  calculated.  These  are  summarised  in 
Table  3. 

The  results  in  Table  3  show  that  the  T91  fuel  cladding  will  pre¬ 
serve  its  integrity  indefinitely  in  an  LBE-cooled  2  GW(th)  ADSR 
when  subjected  to  cyclic  temperature  variations  due  to  pulsed- 
beam  operation.  It  should  be  emphasised  that  the  beam  character¬ 
istics,  in  particular  the  assumed  beam-off  time,  have  been  chosen 
conservatively  so  that  the  thermal  stresses  experienced  by  the 
cladding  are  considerably  more  extreme  than  are  likely  to  be  expe¬ 
rienced  in  practice  if  FFAG  accelerators  capable  of  driving  ADSRs 
are  successfully  developed. 

4.  Conclusions 

Maintaining  the  integrity  of  the  fuel  pin  cladding  during  its  ser¬ 
vice  inside  a  nuclear  reactor  is  of  high  importance  as  the  cladding 
provides  a  barrier  against  radioactive  fission  products  and  prevents 


them  from  leaking  out  into  the  primary  cooling  circuit.  Operating 
an  ADSR  with  a  pulsed  neutron  source  is  one  of  the  possible  sce¬ 
narios  associated  with  the  development  and  deployment  of  FFAG 
accelerator  technology.  The  effect  of  the  temperature  fluctuations 
within  the  cladding  due  to  such  an  operational  regime  was  inves¬ 
tigated  using  the  PTS-ADS  code,  which  was  developed  specifically 
for  this  purpose.  It  has  been  found  that,  even  under  very  conserva¬ 
tive  modelling  assumptions,  the  thermal  fatigue  effects  in  the  clad¬ 
ding  are  very  small  and  therefore  the  number  of  allowable  cycles  to 
fatigue  failure  is  so  high  that  the  mechanical  integrity  of  the  clad¬ 
ding  can  be  assumed  to  be  unaffected  by  the  repetitive  tempera¬ 
ture  variations  caused  by  pulsed-beam  operation  mode. 

5.  Future  work 

The  work  done  in  this  paper  represents  a  first  step  towards  a 
greater  understanding  of  the  reactor-accelerator  coupling  in  ADS¬ 
Rs.  The  PTS-ADS  code  will  be  used  to  perform  a  thermo-mechanical 
analysis  to  determine  the  impact  of  frequent  beam  interruptions 
on  the  mechanical  integrity  of  the  fuel  cladding.  In  contrast  to 
the  situation  in  pulsed-beam  operation,  the  fuel  cladding  is  ex¬ 
pected  to  experience  thermal  shocks  and  therefore  large  stresses 
during  beam  interruptions.  Another  accelerator-induced  transient 
to  be  looked  at  is  beam  overpower  where  a  sudden  increase  in 
the  beam  power  can  also  lead  to  a  thermal  shock  in  the  fuel  clad¬ 
ding.  Additionally,  the  effects  on  cladding  fatigue  damage  of  oper¬ 
ational  and  design  variables,  such  as  irradiation  damage,  high 
temperature  creep  and  choice  of  coolant,  will  be  studied. 

The  pulsed-beam  operation  case  will  be  revisited  to  illustrate 
the  possible  use  of  pulsed  operation  to  measure  the  reactivity  of 
the  core  of  an  ADSR  accurately.  Previous  studies  have  shown  that 
this  is  achievable  by  adopting  methods  like  those  mentioned  in 
(Ryves  and  Scott,  1962)  and  (Baeten  and  Ait  Abderrahim,  2003). 
These  methods  require  a  change  in  the  neutron  flux  to  be  able 
measure  the  reactivity.  It  has  been  suggested  that  this  could  be 
achieved  by  relying  on  accelerator  beam  interruptions.  Pulsed- 
beam  operation,  on  the  other  hand,  has  the  potential  to  provide  on¬ 
line  and  continuous  core  reactivity  measurements  without  the 
need  for  frequent,  potentially  problematic,  beam  interruptions. 
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Appendix  A 

bt=-£yfMT  0 

af 

C1  =4j'Ak/(Ti) 

af 

where  yf.  =  [pf.CPf_ f1  and  df=  rflnf. 
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For  1  <  i  <  nf. 

a>=V-7iyfMT^ 

CXf  1 1 

bt  =  -4  yfMT‘) 


s-xs-'-pyu.w 


_  j_  rff± 2y  U 

n/  "  d}  r„f  Vs 


rnf+fg+a-rnrlJ<f(Tnf) 


where  d}=df  + 1 dg, rn/_,  =  r„f  - 1 df, r„/+,  =  r„f  +  \ds  and  dg=rg-  rf. 

a  —  1  r"/+1^v  h 

a"'+1  “  d*c  r„/+1  yc"/+A 

■ 1  —  —  ^  ^3^cn^+l  rnj+\lhg  +  ^~^nf+^cdnf+l) 

Cn^=idcT^ycnf+MTn^) 

where  ycj  =  [pC(Cpci]_1 ,  dc  =  [rc  -  rg]/nc,  d*c  =  dc+±dg 
r„/+1i  =  rn/+1+idc. 

For  n/+  1  <  i  <  nf+  nc\ 

a.  =  4  rfr-ycMTi) 

dc  n 

&.•  =  -4  ycMTi) 

dr 


and 


1  r,,  i 

Ci  =  -^^ycMTi) 
dc  n 


Qrif+nc  —  \2  ycm+Tic^cy1  nf+ric) 

Clc  1  rif+ric 


brif+ric  ~ 


J__l 
dc  r, 


rif+tic 


ycnf+ric 


rnf+nchw  +  ^  Enf+nc-d<c(Tnf+nc) 
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